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Abstract. We evaluate three typical four-loop non-planar massless propagator diagrams in a Taylor ex- 
pansion in dimensional regularization parameter e = (4 — d)/2 up to transcendentality weight twelve, using 
a recently developed method of one of the present coauthors (R.L.). We observe only multiple zeta values 
in our results. 
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Analytic results for one-scale multiloop Feynman inte- 
grals in a Laurent expansion in e = (4 — d) /2 are expressed 
as linear combinations of transcendental constants with 
rational coefficients. The set of these constants essentially 
depends on the type of Feynman integrals. Probably, the 
simplest type of one-scale Feynman integrals are massless 
propagator integrals depending on one external momen- 
tum. Here the world record is set at four loops — see 
Ref. [I] where all the corresponding master integrals were 
analytically evaluated in an epsilon expansion up to tran- 
scendentality weight seven. 

Practical calculations show that only multiple zeta val- 
ues (MZV) (see, e.g., 0) appear in results. Brown proved 
[3] that convergent scalar massless planar propagator di- 
agrams with the degree of divergence lu = Ah — 2L = —2 
(where h and L are numbers of loops and edges, corre- 
spondingly) up to five loops contain only MZV in their 
epsilon expansions. (In two loops, a proof was earlier pre- 
sented in Ref. 4.) He also proved that for the three non- 
planar diagrams depicted in Fig. [TJ every coefficient in a 
Taylor expansion in e is a rational linear combination of 
MZV and Goncharov's polylogarithms [5] with sixth roots 
of unity as arguments. 

The goal of this brief communication is to study these 
diagrams experimentally. We present results in an epsilon 
expansion up to transcendentality weight twelve. To do 
this we apply the DRA method recently suggested by one 
of the authors (R.L.), Ref. [5J. The method is based on the 
use of dimensional recurrence relations (DRR) [7] and an- 
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Fig. 1. Diagrams considered in Ref. [5] as possible candidates 
for the presence of sixth roots of unity in the e-expansion. 



alytic properties of Feynman integrals as functions of the 
parameter of dimensional regularization, d, and was al- 
ready successfully applied in previous calculations [8ll9lfT0l 
ITHIT2] . To apply this method it is essential to perform an 
integration by parts (IBP) [T3] reduction of integrals that 
participate in dimensional recurrence relations to master 
integrals. To do this, we use the C++ version of the code 
FIRE [Ti] , 

To study analytic properties of solutions of dimen- 
sional recurrence relations, i.e. to reveal the position and 
the order of poles in d in a basic stripe, we used a sector 
decomposition [TMTtjlITT] implemented in the code FIESTA 
[T7lfT8] . To fix remaining constants in the homogenous so- 
lutions of dimensional recurrence relations it was quite 
sufficient for us to use analytic results for the four-loop 
massless propagators master integrals pQ (confirmed nu- 
merically by FIESTA [H]). Finally, after obtaining results 
for master integrals in terms of multiple series we calcu- 
lated resulting coefficients at powers of e numerically with 
a high precision and then applied the PSLQ algorithm 
|20) . We also applied the code HPL [21], a Mathematica 
code dealing with harmonic polylogarithms [35] and MZV. 

For all the three diagrams of Fig. 1, we performed 
evaluation up to transcendentality weight twelve where 
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the basis of transcendental numbers we used includes 48 
constants. Coefficients in our results turn out to be cum- 
bersome, in particular, such are the coefficients at ir 12 in 
subsequent formulae. Therefore, to achieve a successful 
calculation by PSLQ, we were forced to perform numeri- 
cal calculations with a high accuracy. In fact, the accuracy 
of 800 digits was enough to obtain results by PSLQ. How- 
ever, to be on the safe side, we checked our results by 
numerical evaluation with the accuracy of 1500 digits. 

The first diagram of Fig. 1 is a master diagram. This 
is nothing but M45 in Fig. 2 of Ref. [TJ where all the mas- 
ter integrals for four-loop massless propagators are shown. 
Following the method of [6] we needed, first to calculate 
lower master integrals M i, M u , M13, M14, M 2 i, Af 27 , M 35 . 
Eventually, we arrived at the following result which is 
made homogeneously transcendental by pulling out an ap- 
propriate rational function of e and normalizing it at the 
fourth power of the one-loop integral (i.e. M 31 according 
to the notation of [T]): 



non-planar diagrams of Fig. [T] are all planar diagrams. 
However, the result of Brown [3] quoted in the beginning 
of our paper is not applicable to such planar diagrams. 
(Two of them are M 35 and M 36 in Fig. 2 of Ref. [J, with 
uj = 0, rather than —2.) So, a priori, one could admit the 
presence of something additional to MZV in their epsilon 
expansions. 

These are our results for them in the same normal- 
ization where we again managed to reveal homogenous 
transcendentality: 
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are MZV given by 
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The other two diagrams of Fig. 1 are not master inte- 
grals and, consequently, they do not appear in Fig. 2 of 
Ref. [I] . We found their reduction to lower master integrals 
by FIRE and then evaluated resulting master integrals by 
our technique. Let us stress that the master integrals in- 
volved in the IBP reduction of the second and the third 
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We see that only MZV are present in our results. Al- 
though Goncharov's polylogarithms at sixth roots of unity 12 
were allowed to appear according to the analysis of Ref. [3] 
they have not appeared. In fact, such transcendental num.- ^ 
bers do appear in epsilon expansions of some classes of 
massive Feynman integrals 23 . 24 . 251I26J . 

Taking our results into account one could try to prove 
that that there are only MZV in massless propagator di- 
agrams. Unfortunately, the way Brown proceeded [3J in 
the case of a few concrete diagrams can hardly be general- 
ized to other situations because his analysis was based on 
the possibility to perform recursive integrations in Feyn- ^6 
man parameters, as long as the denominator of the inte- 
grand is linear with respect to them, so that some other 
approach is needed. On the other hand, attempts to dis- 17 
cover new types of transcendental numbers, in addition to 
MZV, meet more and more difficulties because we have to 18 
turn to more loops and more terms of e-expansions. Still 
this is certainly possible at the current level of our possi- 
bilities so that we are planning to do this in the nearest 19 
future. 
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